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Abstract 


The  equation 

y"    +   X^p(x)y  =   0 

P(x)    =1  1x1    S   1 

=  a^    1  <  |xl  <  L 
is  investigated.   p(x)  is  taken  as  a  periodic  function  of 
period  2L  and  defined  as  above  in  the  interval  (-L,L). 

The  nature  of  the  characteristic  values  X  for  which  the 

equation  has  solutions  of  period  L  and  2L  is  investigated. 

2 
Stability  charts  are  drawn  showing  for  which  values  of  a  and 

X  all  solutions  are  bounded,  and  for  which  some  solutions  are 
not  bounded. 
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1.   Introduction. 

A  particular   case   of  Hill's   equation  will  be  treated  in  detail,    namely 

y"  +  \   Q(x)y  =  0. 

Here  Q(x)  will  be  taken  as  an  even  periodic  function  of  period  2L  and  de- 
fined by 


Q(x)  =1     0  ^  X  <  1 
=  a^    1  s  X  S  L. 

Evidently  Q(x)  is  a  piecewise  constant  function  with  two  step  discontinui- 
ties in  the  period  interval. 

2.   Preliminary  Remarks. 

We  will  first  quote  some  standard  results  in  the  theory  of  Hill's  equa- 
tion.  Proofs  and  references  may  be  found  in  [l]  and  [2]  . 

When  Q(x)  is  piecewise  continuous,  even  and  has  minimal  period  2L,   we 
can  find  normalized  solutions  with  Initial  values 

y^(0)  =  1   ,   y|(0)  =  0  , 
75(0)  =  0   ,   y^(0)  =  1  . 

The  solutions  of  the  characteristic  equation 

P^  -  2y^(2L)p  +1=0 

are  known  as  characteristic  exponents.   According  to  Floquet's  theorem,  we 


2. 


can  then  find  solutions  f, (x),  fp(x)  such  that 
f^(x  +  2L)  =  P^f^(x) 
f^Cx  +  2L)  =  P2f2(^) 

vhen  P  and  P  are  distinct  characteristic  exponents.   If  p  =  p  ,    i.e. 

P,  =  Pp  =  1  or  P  =  P  =  -  1,  we  can  find  a  solution  y(x)  with  the  property 

y(x  +  2L)  =  Py(x). 

If  P  =  1,  y(x)  has  period  2L;  if  P  =  -  1,  y(x)  has  period  h-L. 

Furthermore,  we  have  a  stability  test .  Solutions  are  stable   if  and 
only  if 

lyi(2L)|  <  1 

or       y^(2L)  =  ±1  aad  y^CSL)  =  y^(2L)  =  0. 

Evidently;,  for  [y  (2L)  |  >  1,  P,  and  P   are  real  and  since  p  p^  =  1  at  least 
one  of  them  is  greater  than  one  in  absolute  value.   The  solution  correspond- 
ing to  that  exponent  evidently  has  increasing  peaJs  values. 

Another  important  statement  is  the  following:   There  exists  a  non-trivial 
periodic  solution  which  is 

1.)   even  and  of  period  2L  if  and  only  if  y^(L)  =  0, 
2.)   odd  and  of  period  2L  if  and  only  if  y^(L)   =  0, 
5.)   even  and  of  period  4l  if  and  only  if  jAl)   =  0, 
h.)      odd  and  of  period  ^L  if  and  only  if  y'(L)  =  0. 

By  means  of  the  Liouville  transfonnation  ^  -'   ve   can  transform  equations 
of  the  fonn 


f2(x)y"(x)  +  f^(x)y'(x)  +  [Xg(x)  +  fQ(x)]y(x)  =  0 


5. 


into  the  form 

Ti"(5)  +  [x  +  Mr(0]Ti(|)  =  0. 

In  order  to  effect  this  transformation  certain  differentiability  conditions 
must  be  imposed  on  the  coefficients  in  the  equation.   In  particular,  to 
apply  the  Licuville  transformation  to  an  equation  of  the  type 

y"  +  X  Q(x)y  =  0 

it  is  necessary  that  Q(x)   possess   a  second  derivative.      In  the  particular 
case  under  discussion  these   conditions   are  not  fulfilled.      A  certain  Oscil- 
lation  Theorem  has  been  proved  by  Haupt    '- -■ ,   which  states   that  the  equation 

where  \|/(|)    is   a  doubly  differentiable,    real  function  of  period  2L,   has 
solutions  of  period  2L  for  an  Infinite   sequence  of  characteristic  values 
X     <  X     <  Xp  <   . .  .    ,    and  solutions   of  period  kl,  for  a  sequence 
X|  <  X'    <  X'    <  . . .    .      These     X's    satisfy  the  Inequalities 

^o  <  ^i  -  ^2  '^  ^1  -  ^2  <  ^5  -  H  <  ^5  -  H  <  •  •  ■ 
and  satisfy  the   relations 

lim       X         =  0   ,    and  '  • 

n 
n  -►   oo 

lim       X'""""  =  0   . 
n 
n  -►  00 

Furthermore,  for  X's  in  the  intervals  (X  ,X'),  (X',X  ),  (XpX'  ),  (X^A  ),  ... 
the  solutions  are  stable  ,  outside  of  these  intervals  unstable  and  at  the 
end  points  they  are  in  general  unstable  ,  except  under  very  special  conditions. 


We  will  investigate  to  what  extent  these  theorems  remain  valid  for 
the  equation  under  dlsciission. 


2 
5.   The  Calculation  of  the  Characteristic  Values  for  a  >  0. 


We  now  turn  to  the  equation 


2 
y"(x)  +  \  y(x)  =0   ,        0  ^  X  <  1 

2  2 
y"(x)  +  X  a  y(x)  =  0  ,   1  S  x  S  L 


and  investigate  the  characteristic  values  under  the  boiindary  conditions 

1.)  y'(0)  =  y'(L)  =  0  (even  solutions  of  period  2L) 

2.)  y(0)  =  y(L)  =  0  (odd  solutions  of  period  2L) 

3.)  y"'(0)  =  y(L)  =  0  (even  solutions  of  period  kh) 

^•)  y(0)  =  y'(L)  =  0  (odd  solutions  of  period  ^L) . 

Corresponding  to  the  first  case,  one  can  show  by  a  direct  calculation  that 

y(x)  =  cos  ^.x   ,   0  S  X  S  1 
=  cos  X  cos  Xa(x  -  1) sin  X  sin  Xa(x  -  1)^,   1  S  x  S  L. 

These  solutions  obviously  satisfy  the  condition  y'(o)  =  0,  and  from 
y'(L)  =  0  we  obtain,  immediately,  the  equation 

a  cos  X  sin  Xa(L  -  1)  +  sin  X  cos  Xa(L  -  1)  =  0. 

An  equivalent  form  is 

tan  X  =  -  a  tan  Xa(L  -  l) . 


In  a  similar  manner  we  can  construct  corresponding  equations  for  the  other 
three  cases  and  can  write  down  the  following  set  of  characteristic  equations: 

1.)   tan  \  =  -  a  tan  Xa(L  -  l) 

2.)   tan  X  =  -  i  tan  Xa(L  -  l) 
'  a       ^ 

5.)   tan  X  =  a  cot  Xa(L  -  l) 
k.)      tan  X  =  -  cot  Xa(L  -  l) . 

We  could  arrive  at  the  characteristic  values  of  X  by  starting  with  the 
characteristic  exponents  defined  by 

P^  -  2y^(2L)p  +1=0, 

and  asking  for  the  values  of  X  for  which 

y^(2L)  =  1  and  y^(2L)  =  -1. 

The   first  of  these  refers  to  the   cases   labeled  1.)    and  2.)    above   and  the 
second  to  the   cases   5.)    and  h,).      A  direct   calculation  shows  that 

y, (2L)   =   cos   2Xa(L  -   l)    cos  2X  -  ^Ca  +  -)    sin  2X  sin  2Xa(L  -   l) 

and  we   can  see  that 

y,(2L)    -   1  =   -   2  cos  X.   co3'^a(L  -   l)  [tan  X  +   a  tan  Xa(L  -   ij] 

X      [tan  X  +  -  tan  Xa(L  -   l)] 


showing  that  the  roots  of 


y^(2L)  -1=0 
correspond  to  the  roots  of  the  chaaracteristlc  equations  of  cases  1.)  and  2.) 
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Similarly, 

y  (2L)  +1=2  cosS^  sinS^a(L  -  l)  [tan  X  -  a  cot  Xa(L  -  l)] 

X   [tan  X  -  i  cot  Xa(L  -  l)] 

showing  the  factorization  corresponding  to  cases  3« )  ^^^   ^•)- 

The  following  method  will  be  used  to  obtain  the  characteristic  values. 
We  treat  case  1.)  as  typical  and  write  the  equation  in  the  form 

Xa(L  -  1)  =  -  tan"  -  tan  X. 
^      '  a 

In  the  evaluation  of  the  inverse  tangent  we  must  take  the  multivalued 
character  of  the  function  into  account.  We  will  graph  both  the  right  side 
and  the  left  side  of  the  equation  and  wherever  the  curves  intersect,  we 
evidently  have  a  characteristic  value. 

Figure  1  was  drawn  for  a  >  1.   For  a  =  1  the  wavy  lines  become  straight 
and  have  slope  -1.   If  we  superimpose  on  this  figure  the  other  three  cases 
we  obtain  fig\ire  2. 

Figure  2  was  drawn  for  a  >  1.   For  a  =  1  the  curves  of  the  type  1.)  and 
2.)  become  identical  straight  lines  of  slope  -1,    and  similarly  for  the 
curves  of  type  5.)  and  k.).     Wherever  the  line  Xa(L  -  l)  =  f(X)  intersects 
these  cixrves,  we  obtain  characteristic  values.   If  -j  X.  I  denote  the  charac- 
teristic )Js  corresponding  to  solutions  of  period  2L  and  |  X'  >    the  ones  corres- 
ponding to  solutions  of  period  k-L,   we  find  that 

0  =  X^  <  X|  i  X^  <  X^  S  Xg  <  X^  ^  Xj^  <  X^  ^  Xi^  <  . . .  . 

Thus  Eaupt's   Oscillation  Theorem  remains  valid. 
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Figiire  1 


I.       2 


Figure  2 


For  a  <  1  the  curves  of  type  1.)  and.  2.)  change  place  and  similarly 
for  types  J.)  and  ^.)' 

We  can  inmedlately  determine  where  multiple  roots  occur.  Evidently 
for  a  =  1  ve  have 


21+1    2i+2 


^"21+1  "  -^21+2 


i  =  0,1,2,.. 


In  general  even  when  a  ^  1,  we  have  multiple  roots  for  a(L  -  l)  rational, 
say  p/q,  where  p  and  q  are  relatively  prime.   Then  we  have  double  roots  at 
the  lattice  points  (kp,kq)  such  that  for 


k(p  +  q)  =  2n  ,       1^_^  =  X^ 


k(p  +  q)  =  2n  +  1  ,    X^^^  =  X^ 


y    n  =  1,2,...  . 


It  also  follows  that  for  a(L  -  l)  =  p/q,  we  need  to  compute  only  the 
first  p  2  roots  of  each  kind  and  all  others  are  equal  to  these  moduio  qn. 
For  a(L  -  1)  =  1  we  can  also  see  that  since 


ten   X^^_^^  =  a  cot  X^^^^ 


*^  ^2i+2  =  i  =°^  ^2i+2 


it  follows  that 


^^^2i+l  ^^^2i+2  =  1 


so  that 

"^21+1  '  "2i+2~  2 


^A.  .-,  +  >^A^  .^-   5     mod  rt. 
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To  obtain  estimates   for  the  roots  we   see  that   for  case  1,) 


1  1 


Xa(L  -  l)   =   -   tan       —  tan  X  «  nrt   -  X. 


From  this  we   see  that 


2n-l 


2n 
1  +   a(L   -   1) 


o  ^    ^    o    •> 


0   S   0   <  1. 


(Here    [x]    represents  the  greatest  integer  less   than  or  equal  to  x. )      Simi- 
larly, 


2n 


2n 
1  +   a(L  -   1) 


^2n-l 


2n  -   1 
1  +   a(L   -   1) 


2+0   2' 


^2n  = 


2n 


_1  +  a(L  -   1)  _ 


2 


To  ohtaln  better  numerical  estimates  in  particular  cases  we  can  use 
Newton's  method. 


k.        The  Stability  Chart  for  a  >  0. 

We  can  make  a  rough  stability  chart  from  Fig.  1.  That  is,  we  will  make 
a  graph  showing  the  relationship  between  the  characteristic  X's  and  a.   These 
curves  subdivide  the  first  quadrant  into  regions  of  stability  and  instability. 
That  is, for  any  X  and  a  from  a  stable  region,  all  solutions  of  the  differential 
are  stable  .   Outside  the  stability  regions  there  is  only  one  bounded  solution. 


u. 


12. 


The  shaded  regions  represent  regions  of  stability.  The  values  of  X 
vhere  the  ctirves  of  type  2.)  and  5-)  intersect  the  abscissa  are  roots  of 
the  equations 

tan  X  =  -  X(L  -   l) 
^^  ^  =  X(L  -   1) 


respectively. 


2 

5.        The   Calciilatlon  of  Characteristic  Values   for  a     <  0- 

2 
To  treat  the   case  where   a     <  0,   we   return  to  the  previous  equations 

and  replace   a  by  ia.      Then  we  have 

y"  +  X^y  =  0   ,  0  g  X  <  1 

y"    -  X^a^y  =  0   ,  1  <  x  g  L  , 


subject  to  the  boiondary  conditions 

1.)  y'(0)  =  y'(L)  =  0 

2.)  y(0)   =  y(L)   =  0 

5.)  y'(0)   =  y(L)   =  0 

^.)  y(o)  =  y'(L)  =  0. 

The  corresponding  characteristic  equations  are 

1.)   tan  X  =  a  tanh  Xa(L  -  l) 

2.)   tan  X  =  -  i  tanh  Xa(L  -  l) 
a 

5. )   tan  X  =  a  coth  Xa(L  -  1) 
h-.)     tan  X  =  -  -  coth  Xa(L  -  1). 
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To  evaluate  the  roots  of  these  equations,  we  write,  as  before,  in 
case  1. ) 


1  1 


Xa(L  -  l)  =  tanh   —  tan  \ 


and  graph  the  right  side  of  the  equation.  We  then  obtain  Flg:ure  k,   where 
all  four  cases  are  included. 


Figure  h 


H^. 


2 

For  a  (L  -  l)  <  1  we  again  have  a  sequence  of  characteristic  roots 

In  the  order  given  by  the  Oscillation  Theorem 


0  =  \  <  \'  <  X'  <  >.,  <  X^  <  X;  <  \,'  <  X^  <  X|,  < 

o    1    2    1    2    5    ^    3"^ 

2 
In  this  case  ve  obtain  no  double  roots.   For  a  (L  -  l)  >  1  a  direct  analysis 

shows  that  ve  have  two  roots  corresponding  to  case  1. )  in  the  interval 

[o,    tan  a]  .   Then  the  sequence  becomes 

0  =  x^"""^  <  x^^^  <  X'  <  X'  <  X,  <  x^  <  x;  <  X,'  <  x^  <  x,,  <  ...  . 

o      o      12125454 

For  a^(L  -  l)  =  1   X^^^  =  X^^^  =  0. 
^  '  00 

For  large  n  we  see  readily  that 

^      -1 
X„     tii  m(  +  tan     a 
2n 

X'      -  jv    nn  +  tan     a 
2n+l  ^ 

^2n+2  ~  ^^  +  ^^'^  -  *^"^  i         ■      ', 
-11 

,  I  n     J.       I     I  -if      _      x  ay 

2n+l 


Figure  5  Is  the  stability  chart  for  this  problem,   where  the  regions  of 
stability  are  shaded. 
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3"     4         2 


Figure  5 


-  i6 


6.  Applications 

The  equation  discussed  in  Sections  2  and  5  has  numerous  physical 

applications.   It  can  be  derived  from  the  problem  of  a  vibrating  string, 

composed  of  two  homogeneous  pieces.  Another  application  is  to  the  propagation 

of  a  wave  in  a  one -dimensional  structure,  where  the  index  of  refraction  is 

an  even  periodic,  piecewise  constant  function. 

Another  application  is  to  resonant  inductance  capacitance  circuits, 

where  the  inductance  or  capacitance  varies  periodically  with  time.   Such 

M         [5] 

equations  were  treated  approximately  by  Carson    and  Minorskj. 
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